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20 Abstract 

21 Partitioning of mineral dust, pollution, smoke, and mixtures using remote sensing 

22 techniques can help improve accuracy of satellite retrievals and assessments of the aerosol 

23 radiative impact on climate. Spectral aerosol optical depth (x) and single scattering albedo (co 0 ) 

24 from Aerosol Robotic Network (AERONET) measurements are used to form absorption [i.e., co 0 

o o 

25 and absorption Angstrom exponent (a a bs)] and size [i.e., extinction Angstrom exponent (a ex t) and 

26 fine mode fraction of x] relationships to infer dominant aerosol types. Using the long-term 

27 AERONET data set (1999-2010), 19 sites are grouped by aerosol type based on known source 

28 regions to: (1) determine the average co„ and a a bs at each site (expanding upon previous work); 

29 (2) perform a sensitivity study on a a b s by varying the spectral to,,; and (3) test the ability of each 

30 absorption and size relationship to distinguish aerosol types. The spectral co 0 averages indicate 

31 slightly more aerosol absorption (i.e., a ().0<6co 0 <0.02 decrease) than in previous work and 

32 optical mixtures of pollution and smoke with dust show stronger absorption than dust alone. 

33 Frequency distributions of a a b s show significant overlap among aerosol type categories and at 

34 least 10% of the a a b s retrievals in each category are below 1.0. Perturbing the spectral co 0 by 

35 ±0.03 induces significant a a b s changes from the unperturbed value by at least -+0.6 for Dust, 

36 -+0.2 for Mixed, and -+0.1 for Urban/Industrial and Biomass Burning. The co 0 440 nm and a ext 440 - 

3 7 870 nm relationship shows the best separation among aerosol type clusters, providing a simple 

38 technique for determining aerosol type from surface- and future space-based instrumentation. 

39 
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40 1.0 Introduction 

41 Particles suspended in the atmosphere are difficult to characterize both temporally and 

42 spatially due to their short lifetime and geographically diverse sources. Aerosol mixtures — 

43 whether dominated by dust, sulfate, carbon, sea salt, or mixtures of these particles — pose a 

44 challenge to satellite and sub-orbital remote sensing techniques when identifying aerosol type 

45 [Jeong and Li 2005; Levy et al., 2007; Kalapureddy et al., 2009; Lee el al., 2010; Kahn et al., 

46 2010; Russell et al., 2010]. Remote sensing techniques can quantify the aerosol particle size 

47 using spectral aerosol optical properties, but inferring aerosol type requires knowledge of the 

48 source regions usually obtained through use of ancillary data sets (e.g., back trajectory models, 

49 satellite product, electron microscopy) to determine emission sources, transport mechanisms, 

50 composition, and morphology. The discrimination of aerosol types increases accuracy of the 

51 assessment of the aerosol radiative impact and therefore is important to climate modeling [Diner- 
Si et al., 1999; Satheesh and Moorthy 2005]. Variations in spectral aerosol absorption magnitudes 

53 can enable partitioning among aerosols from various source regions, fuel types, or combustion 

54 phases. Aerosol absorption together with size can potentially determine dominant aerosol type 

55 from remote sensing and in situ measurements. 

56 Various methods have been proposed using aerosol optical and microphysical properties 

57 to distinguish aerosol types. The magnitude of the aerosol optical depth (AOD, T ex t) and the 

o 

58 spectral dependence of AOD with respect to wavelength (i.e., Angstrom exponent, a ex t) is 

59 commonly used in aerosol remote sensing to infer dominant aerosol types given knowledge of 

60 the source region or typical aerosol transport mechanisms [e.g., Kalapureddy et al., 2009, Boselli 

o 

61 et al., 2012]. Other techniques using the derivative of the Angstrom exponent or spectral 

o 

62 difference of Angstrom exponent wavelength pairs along with aerosol loading and particle 
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63 effective radius may provide further information on particle type with respect to size and growth 

64 of particles [Gobbi et al., 2007; Basart el al., 2009]. Although size varies among particle types, 

65 the spectral absorption also varies. Studies [ Omar et al., 2005, Levy et al., 2007; Mielonen et al., 

66 2009; Lee et al., 2010; Russell et al., 2010] have suggested relationships utilizing the aerosol 

67 absorption and size properties to determine the dominant aerosol type from Aerosol Robotic 

68 Network (AERONET) retrievals [Holben et al., 1998; Dubovik et al., 2000; Dubovik et al., 2002; 

69 Dubovik et al., 2006]. Information content from these relationships varies from generic 

70 identification of major aerosol particle types (e.g., dust, mixed, urban/industrial pollution, and 

71 biomass burning smoke) to specific degrees of absorbing aerosols. Recently, Russell et al. 

o 

72 [2010] have proposed using the absorption Angstrom exponent (AAE, a a b s ), the spectral 

73 absorption aerosol optical depth dependence on wavelength, to further define aerosol type from 

74 AERONET retrievals. For comparison to the Cloud- Aerosol Lidar with Orthogonal Polarization 

75 (CALIOP) instrument, Mielonen et al. [2009] utilized the AERONET single scattering albedo 

76 (co„) difference between 440 and 1020 nm (as suggested by Bergstrom et al. [2002] and 

77 implemented by Derimian et al. [2008]) and a ex t to estimate aerosol type. Further, Lee et al. 

78 [2010] modified this relationship to use only co 0 from 440 nm and the fine mode fraction (r|) of 

79 the AOD at 550 nm to determine the particle size partitioning. Other techniques using spectral 

80 lidar ratios and multiple aerosol optical and microphysical properties retrieved from AERONET 

81 have been implemented to determine aerosol type categories for various applications [Cattrall et 

82 al., 2005; Omar et al., 2005; Qin and Mitchell 2009; Burton et al., 2012]. 

83 In this study, for data available between 1999 and 2010, 19 AERONET sites were 

84 classified by dominant aerosol type [i.e.. Dust, Mixed, Urban/Industrial (U/I), and Biomass 

85 Burning (BB)] based on previous literature. First, aerosol absorption parameters (i.e., co„ and 
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86 a a bs) were analyzed and compared to previous work. Second, sensitivity tests were performed on 

87 the dabs by perturbing co 0 to determine variability within each dominant aerosol type category. 

88 Last, the absorption and size relationships were evaluated and compared to each other based on 

89 the dominant aerosol type categorizations. 

90 2.0 Instrumentation and Method 

91 The Aerosol Robotic Network is a ground-based network of standardized Cimel Sun and 

92 sky scanning radiometers measuring AOD at multiple wavelengths from 340 to 1640 nm and 

93 retrieving other columnar optically effective aerosol properties (e.g., volume size distribution, 

94 complex index of refraction, and single scattering albedo) from sky radiance measurements at 

95 four wavelengths: 440, 675, 870, and 1020 nm [Holben et al., 1998]. The AOD estimated 

96 uncertainty varies spectrally from ±0.01 to ±0.02 with the highest error in the ultraviolet 

97 wavelengths [ Holben et al., 1998; Eck et al., 1999] and calibrated sky radiance measurements 

98 typically have an uncertainty less than 5% [ Holben et al., 1998]. Further descriptions of the 

99 instrumentation, calibration, methodology, data processing, and data quality are described 

100 elsewhere [ Holben et al., 1998; 2006; Eck et al., 1999; 2005; Smirnov et al., 2000]. For all sky 

101 radiance wavelengths (440, 675, 870, and 1020 nm), the co 0 uncertainty is expected to be ±0.03 

102 based on Version 1 almucantar retrieval computations when T44o nm >0.4 [ Holben et al, 1998; Eck 

103 et al., 1999; Dubovik et al., 2002]. When compared to AERONET retrievals, in situ 

104 measurements of co„ were within AERONET uncertainty estimates [Leahy et al., 2007; Johnson 

105 et al., 2009; Muller et al., 2010; Toledano et al., 2011]. 

106 In-depth discussions of the almucantar retrieval products are given by Dubovik and King 

107 [2000] and Dubovik et al. [2000; 2002; 2006] and quality criteria are discussed by Holben et al. 
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108 [2006] . Dubovik et al. [2002] provided averaged almucantar retrieval aerosol optical and 

109 microphysical properties based on aerosol types and source region using AERONET pre-Version 

110 1 data (i.e., data collected and analyzed prior to the release of quality assured Version 1 retrieval 

111 data set in 2003). These results have been used throughout the literature to define aerosol type 

112 based on the aerosol absorption characteristics [ Russell et al., 2010 and references therein]. 

113 Notably, the Version 2 retrievals (i.e., released in 2006) utilized new input data sets (e.g., NCEP 

114 reanalysis, MODIS ecosystem type-based BRDF functions, and geographically and temporally 

115 varying black sky albedo), more dynamic calculations of the surface reflectance than the Version 

116 1 assumption of a green Earth surface reflectance, robust quality checks of the measured sky 

117 radiance inputs, and improved criteria for acceptable sky residual fits [Holben et al., 2006; Leahy 

118 et al., 2007; Eck et al., 2008 and references therein]. For example, in the United Arab Emirates 

119 and Arabian Gulf, Version 2 improvements provided more consistent co 0 magnitudes and spectra 

120 for coarse-mode dust aerosol over two vastly different surfaces (i.e., small island versus bright 

121 desert) with co„ differences of less than 0.01 compared to 0.03 for the Version 1 spheroid 

122 inversion model and with increased absorption at 440 nm, which typically occurs in iron-rich 

123 desert dust, rather than spectrally neutral co 0 from Version 1 retrievals [Eck et al., 2008]. 

124 Additional instrument checks were implemented to assess absorption properties from the 

125 Version 2 almucantar retrievals. To improve the quality of the sky radiance measurements used 

126 for almucantar retrievals, instrument collimator consistency checks were performed to remove 

127 potential artifacts (e.g., induced by spider webs in the tube, excessive dust, or contamination on 

128 the sensor head window due to moisture or dust) in the radiance measurements. The sky 

129 radiance measurements at ±6° azimuth from solar zenith — using the solar aureole and sky gains 

130 for instruments with only Silicon detectors — were required to have a percent difference of less 
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131 than 10% spectrally from 440 to 1020 nm. For Silicon and InGaAs detector instruments (where 
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each detector measures in a different collimator tube), the temperature corrected Silicon and 
InGaAs xioaonm difference (Ax) must be less than Ax iimi t of 0.06/m (where m is the optical air 
mass), which results in a Ax iimi t of 0.03 when m equals 2 and 0.06 for the overhead sun (m=l). 
Collimator consistency checks provide an improved method to further quality assure the Level 
2.0 almucantar retrieval data set. 

Measured aerosol optical depth and computed almucantar retrieval products can be used 

o 

to derive additional aerosol properties. The extinction Angstrom exponent (a ex t) was calculated 
from the spectral dependence of AOD or x ex t with wavelength (k) using equation (1) [Angstrom 
1964]: 


a ext = -dln[xext(k)|/dln(A) ^ 

For a wavelength range between 440 and 870 nm typically using 440, 500, 675, and 870 nm 
AOD — and computed by linear regression of In x versus In k — values near 0 indicate mainly 
coarse mode (radius, r>l pm) aerosol particles, while values near 2 indicate mainly fine or 
accumulation mode (r<l pm) aerosol particles [Holben et al., 1991; Kaufman el ah, 1992; Eck el 
al., 1999; Reid et al., 1999; Holben et al., 2001]. The absorption AOD or x a b s is calculated for 
each wavelength using equation (2): 

Tabs(^) = X ext (X) * [ 1 -COoO)] (2 ) 

[Eck et al., 2010; Russell et al., 2010; Giles et al., 201 la]. Similar to a ext , the spectral 

o 

dependence of x a b s with k on logarithmic scale gives the absorption Angstrom exponent or a a b s in 
equation (3): 


a abs = -dln[x ab s(/4l/dln(A) 


(3) 
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150 Assuming a spectrally constant index of refraction, very small spherical black carbon particles (r 

151 -0.01 pm) can have a AT 1 dependence or a a b s of 1.0 [Bergstrom el al., 2002], while larger 

152 optically effective black carbon particles (e.g., r>0. 1 pm) may have a a b s values below 1.0 for 

153 large cores and up to 1.6 for various shell coatings [Lack and Cappa 2010]. Russell et al. [2010] 

154 analyzed AERONET pre-Version 1 almucantar retrievals from Dubovik et al. [2002] and showed 

155 a a b s values vary between -1.2 and 3.0 for Dust, -0.75 and 1.3 for U/I, and -1.2 and 2.0 for BB. 

156 Eck et al. [2010] analyzed AERONET Version 2 almucantar retrievals and showed sites 

157 dominated by optical mixtures of dust, smoke, and pollution had a a b s values between -1.2 and 

158 1.8 for mixed size particles (i.e., fine mode fraction of the AOD at 675 nm ranged between -0.35 

159 and 0.65). In the present study, the fine mode AOD (if) and coarse mode AOD (t c ) from the 

160 almucantar retrieval — as inferred from the size distribution and refractive indices — were 

161 interpolated to 550 nm using the linear fit of the logarithms of if, x (i.e., Xf+x c ) for 440, 675, and 

162 870 nm wavelengths to calculate the fine mode fraction of the AOD [i.e., r| = Tf/(xf+x c )] at 550 

163 nm (p 550 nm)* 

164 Nineteen AERONET sites were selected for the analysis based on the availability of an 

165 extensive data record (i.e., greater than five data equivalent years of AOD measurements from 

166 1999 to 2010) and the geographic distribution among aerosol source regions (Figure 1). The 

167 sites were designated as one of four commonly used aerosol classifications: Dust, Mixed, 

168 Urban/Industrial (U/I), and Biomass Burning (BB). The classifications were established based 

169 on the source regions and known seasonal changes in aerosol type over these regions (see 

170 references in Table 1). Further, these selected sites should be subject to high aerosol loading 

171 (i.e., T440nm>0.4) to meet the Version 2, Level 2.0 almucantar retrieval sensitivity requirement for 

172 absorption parameters (e.g., a> 0 , "tabs) [Dubovik et al., 2000; Holben et al., 2006]. Sea salt (as well 
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173 as biogenic) aerosols as a dominant particle type category were not considered in this study since 

174 xsoonm is typically less than 0.1 for pure maritime environments [ Smirnov et al., 2002]; however, 

175 for maritime locations affected by aerosol plumes (e.g., Saharan dust transport over Cape Verde 

176 islands), i440nm>0.4 can be satisfied [Smirnov et al., 2009]. Hence, the i440nm>0.4 criterion biases 

177 the data set only to high aerosol loading periods to ensure enough radiometric sensitivity to 

178 compute absorption reliably [Dubovik et al., 2002]. Although Dust, U/I, and BB categories may 

179 represent the dominant aerosol type, episodic aerosol incursions outside of their classification 

180 category likely have occurred at any site during the analysis period (e.g., dust over Shirahama or 

181 Lake_Argyle, biomass burning smoke over GSFC) [ Sana et al., 2003; Qin and Mitchell 2009; 

182 Eck et al., 2003b; O ’Neill et al., 2005]. The Mixed aerosol category encompasses sites primarily 

183 affected by different mixtures of aerosol types (e.g., dust and pollution or dust and biomass 

184 burning smoke mixtures) on a seasonal basis, increasing the probability of at least an optical 

185 mixture state [Derimian et al., 2006; Eck et al., 2010]. Although no explicit seasonal 

186 partitioning is performed, the X440nm>0.4 criterion captures mainly seasonal increases in aerosol 

187 loading at some sites (e.g., GSFC and Mongu) [ Holben et al., 2001]. 

188 3.0 Results 

189 3.1 Retrieved Absorption Properties by Dominant Aerosol Type 

190 Dust particles aggregated with varying combinations of clay, quartz, and hematite exhibit 

191 strong absorption in the blue wavelength region (e.g., 440 nm) with lower absorption in the 

192 visible and near infrared wavelengths (i.e., co 0 increasing with wavelength) \ Sokol ik and Toon 

193 1999]. For fine mode particles (rcl.O pm in the volume size distribution), hygroscopic aerosol 

194 particles (e.g., sulfates) have near neutral spectral dependence and high scattering efficiency 

195 [ Dubovik et al., 2002]. Black carbon (BC) particles have the strongest absorption in the near- 
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196 infrared (co 0 decreasing with X when the sole absorber), while aerosols composed of brown 

197 carbon (BrC) or organic carbon (OC) exhibit stronger absorption in ultraviolet and visible bands 

198 (co 0 increasing with X when the sole absorber) [Eck et al., 2009]. Varying concentrations of BC 

199 particles optically mixed with dust, BrC, and/or OC can produce ambiguous co (l wavelength 

200 dependence (i.e., increasing, decreasing, or constant with X)\ however, the net effect is stronger 

201 absorption across the retrieved spectrum (e.g., 440 to 1020 nm) [Dubovik et al., 2002; Giles et 

202 al, 2011a], 

203 The AERONET Version 2, Level 2.0 absorption properties at each site are presented in 

204 Figure 2 and Table 2 to provide an update to Dubovik et al. [2002] and Russell et al. [2010]. The 

205 spectral co 0 behavior is similar to Dubovik et al. [2002] for most regions. For Solar Village 

206 (Dust), Capo Verde (Dust), GSFC (U/I), Mexico City (U/I), and Mongu (BB), the to,, differences 

207 between Dubovik et al. [2002] and Table 2 (i.e., co 0 Dubovik 2002 - co 0 Table 2 ) showed an overall 

208 average decrease of 0.01 for these sites with the largest decrease of 0.02 spectrally for GSFC and 

209 Capo Verde and smallest decrease ranging from 0 to 0.01 for Mongu. Notably, the co 0 standard 

210 deviations are significantly greater by 0.01 to 0.03 in the present study than Dubovik et al. [2002] 

211 for all five sites. Table 2 differs from Dubovik et al. [2002] due to utilizing different analysis 

212 criteria (e.g., T44o nm >0.4 in Table 2 vs. iio20nm>0.3 and a ex t<0.6 for desert dust in Dubovik et al. 

213 [2002]), implementing improved surface characterization and inversion quality checks in 

214 Version 2 (as discussed in Section 2), and utilizing a larger data set (e.g., the number of ro 0 

215 retrievals at GSFC is four times larger than Dubovik et al. [2002]). For 0 O 44 Onm as a function of 

216 i 440 nm, the R 2 values — calculated based on a second order fit — ranged from 0.0 to 0.16 for each 

217 site, indicating weak correlation and only up to 16% of a> O 440 nm variation was explained by T 44 o nm - 

218 Table 2 shows that the Dust category has the least variability among sites likely due to the 
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219 similar mineral composition, while the BB category has the largest variability likely due to 

220 various fuel types and fuel combustion phases resulting from different relative BC emissions 

221 [Eck et al., 2003b]. The Mixed category (0.33<r) 5 5o n m<0.66) co 0 average shows strong spectral 

222 absorption and dust-like co 0 spectra with stronger absorption at 440 nm due to significant dust 

223 contribution to the optical mixture. Sokolik and Toon [1998] showed that varying hematite 

224 amounts in dust can lead to increased absorption spectrally from the blue to near-infrared 

225 wavelength region. Using a ex t<0.2 to designate “pure dust” as suggested by Kim el al. [201 1], 

226 the overall “pure dust” average of co 0 for all Dust category sites is 0.91, 0.97, 0.97, 0.97 for the 

227 440, 675, 870, and 1020 nm wavelengths, respectively. These “pure dust” co 0 values are lower 

228 by up to 0.02, spectrally, than those reported by Dubovik et al. [2002] for Dust sites and are 

229 lower by up to 0.01 for oo 0 at 550 nm (logarithmically interpolated between 440 nm and 675 nm) 

230 compared to similar sites analyzed by Kim et al. [201 1]. Table 2 shows the Dust site averages 

231 are lower than “pure dust,” indicating possible incursions by other aerosols (e.g., biomass 

232 burning smoke). Eck et al. [2010] and Giles et al. [201 la] also showed increasing absorption 

233 with wavelength for decreasing a a b s , indicating an optical mixture and possibly aggregation of 

234 dust and carbonaceous particles at Kanpur, India. We interpolated r| to 550 nm using the linear 

235 fit of the logarithms of if, x (i.e., if + x c ) and the 440, 675, and 870 nm wavelengths similar to Lee 

236 et al. [2010]. In Figure 3a, the Mixed category for the coarse mode particles (r| 55 onm: 0.0-0.33) 

237 resembles dust co 0 spectra as shown in Figure 2a. In Figure 3c, for fine mode particles (r) 55 o nm : 

238 0.66-1.0), the co 0 magnitudes and variability are similar to U/I or BB particle types categories but 

239 with less co 0 spectral dependence possibly due to varying amounts of BC, BrC, and OC 

240 [Derimian et al., 2006; Eck et al., 2009; 2010]. The average co„ for a a bs binned between 1.5 and 

241 2.0 shown by Giles et al. [201 la] at Kanpur closely resembles the absorption magnitude and 
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242 spectral shape of mixed aerosol types for the Mixed category (0.33<r|550nm<0.66) in Table 2 as 

243 well as Figure 2b and Figure 3b suggesting various mixtures of aerosol particles contributing to 

244 the absorption. 

245 The co 0 and x ex t are used to derive the x a bs from AERONET data. x a b s and a a b s were 

246 averaged for each site based on dominant particle type category in Figure 4 and Table 2. For the 

247 five sites (Solar Village, Capo Verde, GSFC, Mexico City, and Mongu), a comparison of 

248 average a a b s values in Table 2 with Russell el al. [2010] for the 440-870 nm range shows the 

249 largest difference in a abs (i.e., a abs Ru sseii 2010 - a abs Table 2 ) at GSFC (-0.25) and Capo Verde (+1.2). 

250 For the other three sites, the a a b s averages in Figure 4 are comparable to those reported by Russell 

251 et al. [2010] and Giles et al. [201 la]. In Figure 5, the Mixed category was further stratified by 

252 the r| 550 nm as in Figure 3. The coarse particle range (r) 550 n m : 0.0-0.33) shows similar a a b s (1.7-2. 3) 

253 as the Dust category (which is expected for dust dominated cases) and the fine particle range 

254 0l55Onm- 0.66-1.0) shows an a a b S (0.8- 1.5) similar to BB and U/I categories. The mixed size 

255 particle range (r| 550 nm: 0.33-0.66) is nearly identical to the Mixed category a a b s (1.2- 1.7) in 

256 Figure 4b and similar to values reported by Eck et al. [2010]. As shown by Bergstrom et al. 

257 [2007] and Russell et al. [2010], the a a b s may vary significantly when considering the aerosol 

258 particle size between fine and coarse modes; however, when considering U/I and BB aerosols 

259 within the fine particle range, significant overlap results in a a b s . The sensitivity of a a b s with 

260 respect to input parameters will be investigated in the next section. 

o 

261 3.2 Absorption Angstrom Exponent Sensitivity Study 

262 The retrieved values of a a b s have a normal distribution (Figure 6) when calculating a a b s 

263 using three wavelengths (440-675-870 nm) for each dominant aerosol type. Russell et al. [2010] 

264 showed that the average a a b s values generally decreased with increasing spectral range possibly 
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265 due in part to the crude surface reflectance assumption made in early AERONET analysis (as 

266 discussed in Section 3.1), while Gyawali et al. [2012] showed an increasing a a b s values with 

267 increasing spectral range for clean and polluted days during winter in Reno, Nevada. However, 

268 increasing or decreasing trends of a a b s depend on the wavelength interval [Lack and Cappa 

269 2010]. These a a b s values computed from 440-675-870 nm wavelength range have large 

270 variability with standard deviations ranging from ±0.3 to ±0.6 (1.76±0.58 for Dust; 1.53±0.44 for 

271 Mixed; 1.21±0.37 for U/I; 1.35±0.35 for BB). Individual a a b s retrieval calculations [a a b s (Dust): 

272 -0-4; a a b s (Mixed): -0-3.5; a a b s (U/I): -0-2; a a b s (BB): -0-2.5] are within the range of all 

273 dominant particle types; therefore, a a b s should not be used alone to determine aerosol types 

274 without the use of other information (e.g., aerosol size). Further, Figure 6 shows a significant 

275 percentage of a a b s below 1.0, which is the black carbon limit for very small particles [Bergstrom 

276 et al., 2002]. However, Lack and Cappa [2010] suggested a a b s (from 380-750 nm) values for 

277 larger optically effective BC particles may exist between —0.2 and 1.6 depending on the BC 

278 coating material. Nonetheless, the U/I category has over 22% of the a a b s retrievals below 1.0, 

279 while the other categories have -10% of the data below a a b s of 1.0 also possibly due to the 

280 uncertainty of the retrieved co 0 . 

281 A sensitivity study of a a b s was performed to test the response of a a b s by varying co 0 for 

282 each wavelength (i.e., 440, 675, and 870 nm) and holding x ext constant (i.e., AOD cannot be 

283 perturbed since it is used in the calculation of a> 0 retrieval). First, co 0 was varied by ±0.03, the 

284 current AERONET co 0 uncertainty estimate [Dubovik et al., 2002]. Additional tests included 

285 varying co 0 by ±0.01, ±0.02, and ±0.04 to show the variability of a a b s with various degrees of co 0 

286 uncertainty. Different spectral co 0 inputs schemes were implemented to determine the a a b s 

287 response by varying co 0 equally across all wavelengths, by perturbing co 0 at only one end point in 
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288 the 440-675-870 nm wavelength set (i.e., 440 nm or 870 nm), and by perturbing co 0 at 440 nm or 

289 870 nm in the 440-870 nm wavelength pair (i.e., excluding 675 nm). Positive co (l perturbation 

290 may approach values of 1.0 (i.e., absolute scattering) and can produce large positive or negative 

291 a a bs due to very low x a b s . To prevent such cases, the 0 O magnitude was limited to less than 0.995 

292 for positive co 0 perturbations for all wavelengths resulting in a reduced data subset. 

293 Table 3 shows the sensitivity of a a b s to perturbations in to,,. The perturbation of ±0.03 co 0 

294 (the current estimated uncertainty) changed a a b s by at least ~±0.6 for Dust, ~±0.2 for Mixed, and 

295 ~±0.1 for U/I and BB. The perturbations of co 0 by ±0.02 showed ~0.1 smaller corresponding 

296 change in a a b s with respect to ±0.03 co 0 perturbations for Dust and less than 0.05-0.10 for the 

297 other categories. Perturbations of co 0 by ±0.04 showed large deviations from the unperturbed 

298 data set, indicating much greater uncertainty for a a b s with increasing co 0 uncertainty. The 

299 simulated overestimation of spectral ® 0 for U/I and BB (i.e., 5 cq o =- 0.03) showed a higher 8a a b s 

300 suggesting a possibility that the unperturbed a a b s is underestimated and may possibly, at least 

301 partly, explain a a b s below 1.0 in these categories. However, Lack and Cappa [2010] determined 

302 that the large a a b s variation (-0.2 and 1.3 for the 380-750 nm wavelength range) for BC particles 

303 with coatings are possible and a a b s values less than 1.0 may occur with larger BC particles (e.g., 

304 r core >0.1 pm and r she n>0.25 pm). Gyawali et al. [2012] showed laboratory measurements of 

305 kerosene soot particles have a a b s values of ~0.8 for the 355-1020 nm range and in situ 

306 measurement values of a a b s measured during the Reno, Nevada, winter period varied for clean 

307 days (PM2.5<40 pg/m 3 ) between -1.0 and 1.4 and for polluted days (PM 2 .5>40 pg/m 3 ) between 

308 0.9 and 1.2 for the 405-870 nm wavelength range. Although these model simulations and 

309 laboratory and in situ measurements suggest a a b s values may occur below 1.0, AERONET 

310 remotely sensed values of a a b s have not yet been compared to coincident column-effective in situ 
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311 measurements (e.g., measured by aircraft) but this analysis will be addressed in future work. In 

312 the present analysis, the simulated underestimation of spectral co 0 (i.e., 8co 0 =+0.03) for Dust and 

313 Mixed indicates possible underestimation of the unperturbed a a b s , which could also result in a a b s 

314 below 1.0. Table 3 also shows the Dust and Mixed categories change in the same direction as 

315 the co 0 perturbation possibly due to weak spectral dependence of x ex t, while U/I and BB 

316 categories have the opposite response. Two additional tests were conducted by perturbing to,, 

317 using the wavelength pair (440-870 nm) and only varying the end points of the 440-675-870 nm 

318 set and the differences between unperturbed a a b s averages were minimal (not shown). However, 

319 perturbing one co 0 end point for either the 440-870 nm wavelength pair (not shown) or the 440- 

320 675-870 nm set (Table 3) produced very large deviations in a a b s by up to ~1.2 for Dust, -0.7 for 

321 Mixed, -1.0 for U/I, and -0.6 for BB. The perturbation of end points simulates atypical behavior 

322 of the instrument while deployed in the field (e.g., anomalous filter degradation) showing 

323 potential issues in using real-time data products unless further screening is implemented, such as 

324 the instrument collimator consistency checks (stated in Section 2.0), which may be utilized to 

325 help remove co„ artifacts (i.e., collimator or sensor head window obstructions) and improve the 

326 reliability of a a b s retrievals. These sensitivity tests quantified the effect of the reduction of co n 

327 uncertainty on improving estimates of a a b s . 

328 3.3 Cluster Analysis by Dominant Aerosol Type 

329 Knowledge of aerosol particle spectral absorption provides insight to determine aerosol 

330 particle dominance of dust, carbonaceous matter, or hygroscopic aerosols (e.g., sulfates, nitrates, 

331 or sea salt). While the co-albedo (or 1 -co 0 ) indicates the magnitude of absorption and a a b s 

332 provides some indication of the dominance of carbonaceous particles (e.g., BC, BrC, and OC) or 

333 iron oxides in dust, these parameters alone cannot fully describe the aerosol particle type. 
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334 Recent studies have suggested applying an aerosol particle size parameter (e.g., a ex t or r|) to 

335 separate larger dust particles from other aerosol types and mixtures [Lee et al., 2010; Russell el 

336 al., 2010; Giles et al., 2010; 2011a; 2011b]. In this section, several years of AERONET 

337 retrievals of 0 O 44 Onm, a a bs 440 - 870 nm, a ex t 440 - 870 nm, and r) 550 nm (with wavelength subscripts removed 

338 hereafter) were analyzed for each dominant aerosol type category using a density based 

339 clustering utilizing the Voronoi tessellation [Voronoi 1908; Ishimoto et al., 2010] to determine 

340 the relative concentration of points (density = 1/polygon area) for each absorption and size 

341 relationship. In these density plots (e.g.. Figures 7-10), the high density represents the primary 

342 mode for the dominant aerosol particle type category. Various clustering techniques were 

343 attempted previously to categorize dominant aerosol particle type at AERONET sites [Cattrall et 

344 al., 2005; Omar et al., 2005; Levy et al., 2007; Qin and Mitchell 2009; Russell et al., 2010, 

345 Boselli et al., 2012]. For each absorption and size relationship and aerosol type category in this 

346 study (Figure 11), dominant aerosol particle clusters were computed using averages weighted by 

347 density magnitudes normalized to a 64-level scale (corresponding to a 64-bit color scale). 

348 Although weighting reduces the bias introduced by outliers affecting the normal average, 

349 additional thresholds were applied to the aerosol size parameters. To further define weighted 

350 cluster averages, the a ext cluster averages utilized a 0.8 threshold, where >0.8 indicates mainly 

351 small sub-micron radius particles and <0.8 is mainly large super-micron radius particles (where 

352 a ex t440-870nm = 0 . 8 is approximately equivalent to rpoonm=0.5 as shown for example by Eck et al. 

353 [2005; 2010]). In addition, the q cluster averages were defined using thresholds of 0.0 to <0.33 

354 (for coarse mode dominated particles), 0.33 to <0.66 (for mixed size), and >0.66 (for fine mode 

355 dominated particles). For the BB category (Figure 1 lc and Figure 1 Id), cluster separation was 
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356 imposed to calculate two additional clusters using a co 0 440 nm threshold of 0.90 based on the 

357 density cluster analyses shown in Figure 9 and Figure 10. 

358 The relationships of aerosol absorption and size are analyzed with respect to the 

359 dominant aerosol type category. For Figures 7-10, the primary density clusters are clear 

360 (denoted by orange and red regions representing relative value levels of ~45 to 64). For 

361 example, the Dust category shows a cluster in the region with a ex t of -0.2-0. 3 and p of -0.2-0. 3, 

362 indicative of domination by coarse mode particles. To provide a better assessment of the 

363 clusters, the weighted cluster average and one standard deviation was calculated for each 

364 parameter shown in Figure 11. In Figure 1 la and Figure 1 lb, the primary Dust clusters show 

365 variation of the a a b s mainly between 1.5 and 2.3, which are slightly lower values than reported by 

366 Russell et al. [2010]. In Figure 1 lc and Figure 1 Id, co G also varies significantly in the primary 

367 Dust cluster from 0.89-0.93, possibly due to variation in mineral composition of dust [Sokolik 

368 and Toon 1999]. For a ex t (Figure 11a and Figure 1 lc), the Mixed aerosol type category has two 

369 primary density clusters (1) “Mixed-Large Particle” cluster for mainly super-micron particles 

370 (centered at a ext ~0.4) and (2) “Mixed-Small Particle” for mainly submicron particles (centered at 

371 a ext ~1.25). In comparison to the Dust cluster, the Mixed-Large Particle cluster tends to have a 

372 slightly smaller contribution to larger particles in the 0.3-0. 6 a ext range, while p relationships 

373 (Figure 1 lb and Figure 1 Id) show the Mixed-Large Particle cluster for coarse particles is nearly 

374 identical to the Dust cluster. The Mixed category for mixed sizes (0.33<p55o n m<0.66) does not 

375 show high cluster density due to varying sizes and contributions of the aerosol particles 

376 containing dust with pollution or biomass burning smoke with strongly varying absorption [Eck 

377 et al ., 2010]. In Figure 11, the Mixed-Small Particle clusters (a ex t~1.0 to 1.5; also p~0.8-0.95) 

378 show significant variability likely due to variation in carbonaceous particle contribution 
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379 (primarily BC but also OC) with a a bs between ~1.3 and 1.7 similar to those observed in Kanpur 

380 for fine mode dominated cases by Giles et al. [2011a]. As indicated by Russell et al. [2010] and 

381 shown in Figure 6, the U/I and BB category types for the a a bs and a ex t relationship tend to overlap 

382 each other. For primary density clusters in these two categories, the a a b s vary from -1.1 to 1.8. 

383 Until the co 0 uncertainty is known and constrained further (given the sensitivity results of Section 

384 3.2), the usefulness of a a b s to determine various carbonaceous aerosol particles is doubtful except 

385 in separating cases dominated by BC from cases dominated by BrC or OC. A “region” of higher 

386 a a bs values from the density cluster analysis for BB (Figure 7d and Figure 8d) likely indicates 

387 aged smoke from primarily smoldering combustion containing higher concentrations of BrC or 

388 OC and relatively low BC [Eck et al., 2009; MoosmUller et al., 2009; 201 1], especially above an 

389 a a b s of 1.6 for fine mode particles [Lack and Cappa 2010] and also supported by Figure 9d and 

390 Figure lOd with co 0 above 0.90. For example, according to Eck et al. [2009] and Arola et cd. 

391 [201 1], significant absorbing OC concentrations and high OC/BC ratios likely occurred at the 

392 Bonanza Creek site where Table 2 shows the spectral to,, average is -0.95 and averages of a a b s , 

393 a ext , and r| 550 nm are 1.8, 1.5, and 0.96, respectively. However, the co„ relationships (Figure 11c 

394 and Figure 1 Id) show more cluster separation than a a b s relationships (Figure 11a and Figure 

395 1 lb). In Figure 1 Id, the primary U/I cluster is centered above 0.95, while the main BB cluster is 

396 centered on -0.89 with two BB sub-clusters centered on -0.93 and -0.87 m 0 (calculated by using 

397 0 O threshold of 0.90); however, the BB clusters overlap with the Mixed-Small Particle cluster. 

398 The M O 440 nm and a ext 440-870 nm relationship (Figure 1 lc) shows distinct high density clusters in all 

399 categories (i.e., between Dust and Mixed-Large Particle, and among Mixed-Small Particle-U/I- 

400 BB, and between U/I-BB clusters), while minimal overlap occurs with the U/I and the less 

401 absorbing (m o >0.90) BB sub-cluster. Nonetheless, the analysis has shown that the M O 440 nm and 
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402 a ex t 440-870 nm relationship demonstrates that the dominant particle type may be ascertained simply 

403 from commonly measured or retrieved absorption and size parameters. 

404 4.0 Conclusions 

405 In this paper, absorption properties (i.e., single scattering albedo and absorption 

o 

406 Angstrom exponent) were averaged for 19 AERO NET sites to show correspondence to 

o 

407 representative aerosol source regions. Sensitivity tests on absorption Angstrom exponent were 

408 performed by varying the single scattering albedo within plausible constraints based on 

409 uncertainty estimates. Lastly, the absorption and size relationships were evaluated and compared 

410 to each other based on the dominant aerosol type categorizations. 

411 (1) A summary of aerosol absorption parameters from the AERONET Version 2, Level 

412 2.0 almucantar retrievals was presented to expand upon previous work using pre-Version 1 

413 retrievals. A comparison of five sites common to Dubovik et al. [2002] showed a 0.01 average 

414 spectral (from 440 to 1020 nm) decrease in single scattering albedo (oj 0 ) with the largest 

415 decreases spectrally of 0.02 at Capo Verde and GSFC AERONET sites. The average absorption 

o 

416 Angstrom exponent (a a b S 440-870 nm) computed from Version 2 retrievals was 1.2 lower for Capo 

417 Verde and 0.25 higher for GSFC than reported by Russell et al. [2010] computed from pre- 

418 Version 1 retrievals. Aerosol mixtures exhibited stronger spectral absorption (i.e., lower co 0 ) and 

419 increased dominance of absorbing carbonaceous particles (i.e., lower a a b S 440-870 nm) than for dust 

420 alone, possibly due to an optical mixture state (e.g., dust and smoke or dust and pollution) or the 

421 aggregation of dust and carbonaceous particles. 

422 (2) The a a b S 44 o- 870 nm calculated from AERONET data ranged from ~0 to 3.5 among 

423 dominant aerosol type categories. Frequency distributions of a a b S 44 o -870 nm exhibited significant 
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424 overlap among aerosol types, while the Urban/Industrial and Biomass Burning distributions were 

425 nearly identical for a a bs 440-870 nm values above 1.0. Further, frequency distributions showed 

426 approximately 10% of the a a bs retrievals had values below 1.0 for most aerosol categories but as 

427 high as 22% for the Urban/Industrial category. 

428 (3) A sensitivity study perturbing the co Q by the current AERONET uncertainty (±0.03) 

429 showed a a b s changes by at least -+0.6 for Dust, -+0.2 for Mixed, and -+0.1 for Urban/Industrial 

430 and Biomass Burning. The sensitivity study quantified the improvement in estimates of a a b s 

431 resulting from reducing the co 0 uncertainty. Variations within the uncertainty of co 0 retrievals 

432 may explain some of the observed a a b s values below 1.0 in AERONET data although in situ 

433 measurements suggest that some of these a a b s values may be real depending on the aerosol 

434 particle composition and size. 

435 (4) Absorption and size relationships were examined using density cluster analysis for 

436 each dominant aerosol particle type. The co 0 440 nm vs. a ex t440-870 nm relationship showed at least five 

437 distinct aerosol type clusters [Dust, Mixed-Large Particle, Mixed-Small Particle, 

438 Urban/Industrial, and Biomass Burning (with two sub-clusters)], while the a a b S 440 - 870 nm vs. a ext 440 - 

439 870 nm relationship had fewer distinct clusters due to less definition for mainly small aerosol 

440 particles (a ext4 40-870nm> 1-5). 

441 We showed the co 0 440nm and a ext 440-870nm relationship provided a better clustering 

442 relationship and it may be applied to measurements of aerosol absorption and size properties 

443 derived from surface- and potentially future space-based platforms. From co 0 440 nm and a ext 440-870 

444 nm clusters, at least, major dominant aerosols types and some mixtures can be identified using 

445 common aerosol absorption and size parameters without prior knowledge of aerosol transport or 
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446 source regions. Alternatively, when co 0 is not available but a a b S 440-870nm is (e.g., using spectral 

447 absorption coefficients measured in situ to calculate a a b s ), a a b S 44c-870nm vs. a ex t440-870nm may also 

448 provide a reasonable aerosol type classification. A combination of M O 440nm vs. a ex t440-870nm and 

449 a a bs440-870nm vs. a ext 440-870nm relationships could provide a more detailed classification of aerosol 

450 composition which will require further investigation. 
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686 Table 1 . Previous studies identifying regional aerosol sources affecting AERONET sites 


Aerosol Type 
Source Regions 

Affected AERONET Sites 

Selected References 

Most Regions with 
Various Types 

Most Sites 

Holben et cd. [2001] 

Dust - African 
Dust - Asian 

Banizoumbou, Capo Verde, 
Dakar, Ouagadougou, 
XiangHe, Shirahama 

Tanre et al. [2001]; Reid et al. [2003] 
Eck et al. [2005] 

Smoke -Amazonia 
Smoke - Australian 
Smoke - Boreal 
Smoke - African 

Abracos Hill, Alta Floresta 
Lake_Argyle 
Bonanza Creek 
Mongu 

Eck et al. [2003b]; Schafer et al. [2008] 
Mitchell et al. [2006] 

Eck et al. [2009] 

Eck et al. [2003a; 2003b] 

Pollution - Europe 

Ispra 

Melin and Ziborcli [2005] 

Mixed - Asia 

Mixed - India 
Mixed - Africa 

XiangHe, SEDE_BOKER 

Kanpur 

Ilorin 

Derimian et cd. [2006]; Eck et cd. [2010]; 
Yang et cd. [2009] 

Dey et cd. [2004]; Singh et al. [2004]; 
Prasacl et al. [2007]; Giles et al. [2011a] 
Eck et cd. [2010] 
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687 Table 2. Average aerosol absorption and size properties by aerosol type category from AERONET Version 2 almucantar retrievals 3 


Site 

Date Range 

®o 

440/675/870/1020 nm 

C^abs440-870nm 

C^ext440-870nm 

h550nm 

N 


Dust 


Banizoumbou 

1999-2010 

0.91/0.95/0.96/0.96 

0.04/0.04/0.04/0.04 

1. 7+0.6 

0.3+0. 2 

0.28+0.20 

2901 

Capo_ Verde 

1999-2010 

0.91/0.96/0.97/0.97 

0.03/0.03/0.03/0.03 

2.0+0. 6 

0.2+0. 2 

0.24+0.16 

1202 

Dakar 

2000-2010 

0.89/0.95/0.96/0.96 

0.03/0.04/0.04/0.03 

1 .9+0.6 

0.3+0. 2 

0.28+0.23 

2250 

Ouagadougou 

1999-2007 

0.90/0.94/0.95/0.95 

0.04/0.04/0.04/0.03 

1. 6+0.5 

0.3+0. 2 

0.30+0.21 

1497 

Solar_Village 

1999-2010 

0.91/0.95/0.96/0.96 

0.02/0.02/0.02/0.02 

1. 8+0.6 

0.3+0. 3 

0.28+0.25 

3029 


Mixed (for 0.33<r|55o n m<0.66) 


T1 • 1QQQ onno 0.86/0.90/0.92/0.92 

Ilorin 1999-2009 0 05/0 05/0 04/0 04 1. 6±0.4 0.7+0. 2 0.47+0.23 798 


Kanpur 

2001-2010 

0.87/0.90/0.92/0.93 

0.03/0.03/0.03/0.03 

1. 4+0.4 

0.7+0. 2 

0.48+0.22 

963 

SEDE_BOKER 

1999-2010 

0.91/0.93/0.93/0.94 

0.02/0.02/0.03/0.03 

1. 2+0.5 

0.7+0. 2 

0.48+0.20 

170 

XiangHe 

2001, 2004-2010 

0.88/0.92/0.93/0.93 

0.03/0.03/0.03/0.03 

1. 8+0.4 

0.8+0. 2 

0.53+0.22 

446 

688 a Aerosol optical depth (AOD) at 440 nm is greater than 0.4 for Version 2, Level 2.0 almucantar retrievals. The spectral single 


o 

689 scattering albedo (co 0 ) averages are listed first followed by their standard deviations. The absorption and extinction Angstrom 

690 exponents (a a bs and a ex t) and are computed using the 440-675-870 nm wavelength interval. The fine mode fraction of the AOD 

691 (r| 550 nm) is interpolated to 550 nm as discussed in Section 2. 

692 
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693 


Table 2. (Continued) 


Site 

Date Range 

Wo 

440/675/870/1020 nm 

C^abs440-870nm 

C^ext440-870nm 

d550nra 

N 



Urban/Industrial 





GSFC 

1999-2010 

0.96/0.95/0.94/0.93 

0.02/0.02/0.03/0.03 

1.1 ±0.2 

1. 8+0.2 

0.94+0.20 

882 

Ispra 

1999-2010 

0.93/0.93/0.92/0.91 

0.03/0.04/0.04/0.04 

1. 4+0.4 

1. 6+0.2 

0.92+0.24 

583 

Mexico_City 

1999-2010 

0.89/0.88/0.86/0.85 

0.04/0.04/0.05/0.06 

1. 3+0.3 

1. 6+0.2 

0.87+0.18 

540 

Moldova 

1999-2010 

0.93/0.92/0.90/0.89 

0.03/0.04/0.05/0.05 

1 .2+0.3 

1. 6+0.3 

0.87+0.28 

558 

Shirahama 

2000-2010 

0.94/0.93/0.92/0.92 

0.03/0.03/0.04/0.05 

1.1 +0.5 

1.3+0. 3 

0.81+0.35 

726 



Biomass Burning 





Abracos_Hill 

1999-2005 

0.93/0.91/0.90/0.88 

0.02/0.03/0.04/0.05 

1. 3+0.4 

2.0+0. 1 

0.95+0.14 

342 

Alta_Floresta 

1999-2010 

0.93/0.92/0.90/0.89 

0.02/0.03/0.04/0.05 

1. 5+0.4 

1. 9+0.2 

0.92+0.18 

593 

Bonanaza Creek 

1999-2005, 2008- 
2010 

0.95/0.96/0.96/0.95 

0.03/0.03/0.04/0.04 

1.8+0. 5 

1.5+0. 2 

0.96+0.22 

144 

Lake_Argyle 

2002-2006, 

2009-2010 

0.85/0.83/0.82/0.81 

0.04/0.05/0.06/0.07 

1. 4+0.3 

1.5+0. 4 

0.79+0.36 

176 

Mongu 

1999-2007, 2009 

0.87/0.83/0.80/0.77 

0.03/0.04/0.04/0.05 

1.2+0. 2 

1.9+0. 1 

0.92+0.10 

1411 


694 
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695 Table 3. Sensitivity of the absorption Angstrom exponent (a a b s ) to perturbations of single 

696 scattering albedo (co 0 ) for each dominant aerosol particle type 


Type 

A (nm) 

„ a 

Cl a bs 



5a abs b 


N 

5co 0 =0.0 

OCOo 

All x(A) 

U40nm 

^870nm 



1.76+0.58 

-0.01 

-0.27 



10879 



1.67+0.52 

+0.01 c 

+0.40 



9807 



1.76+0.58 

-0.02 

-0.45 



10879 

Dust 

440-675-870 

1.49+0.42 

+0.02 c 

+0.67 



7290 

1.76+0.58 

-0.03 d 

-0.57 

+0.47 

-0.90 

10879 



1.33+0.38 

+0.03 cd 

+0.79 

-0.54 

+1.16 

4898 



1.76+0.58 

-0.04 

-0.67 



10879 



1.23+0.36 

+0.04° 

+0.85 



3342 



1.53+0.44 

-0.01 

-0.09 



7199 



1.52+0.42 

+0.01° 

+0.13 



7051 



1.53+0.44 

-0.02 

-0.16 



7199 

Mixed 

440-675-870 

1.47+0.38 

1.53+0.44 

+0.02° 

-0.03 d 

+0.23 

-0.21 

+0.40 

-0.53 

6623 

7199 



1.43+0.35 

+0.03 cd 

+0.30 

-0.51 

+0.71 

6060 



1.53+0.44 

-0.04 

-0.25 



7199 



1.40+0.33 

+0.04° 

+0.35 



5479 



1.21+0.37 

-0.01 

+0.05 



3289 



1.20+0.36 

+0.01 c 

-0.10 



3174 



1.21+0.37 

-0.02 

+0.09 



3289 

Urban/ 

440-675-870 

1.19+0.35 

+0.02 c 

-0.21 



2874 

Industrial 

1.21+0.37 

-0.03 d 

+0.12 

+0.74 

-0.52 

3289 



1.18+0.34 

+0.03 cd 

-0.31 

-1.02 

+0.58 

2428 



1.21+0.37 

-0.04 

+0.14 



3289 



1.18+0.34 

+0.04° 

-0.40 



2027 



1.35+0.35 

-0.01 

+0.03 



2666 



1.34+0.34 

+0.01 c 

-0.04 



2639 



1.35+0.35 

-0.02 

+0.06 



2666 

Biomass 

440-675-870 

1.33+0.32 

+0.02 c 

-0.10 



2598 

Burning 

1.35+0.35 

-0.03 d 

+0.08 

+0.45 

-0.31 

2666 



1.32+0.31 

+0.03 cd 

-0.19 

-0.62 

+0.35 

2512 



1.35+0.35 

-0.04 

+0.11 



2666 



1.31+0.29 

+0.04° 

-0.29 



2421 


697 a indicates the unperturbed a a bs average is recalculated based on available t%- 

698 b indicates wavelength(s) used in perturbation of co 0 . 

699 c indicates positive perturbation of co 0 must be less than 0.995 for any wavelength. 

700 d indicates these criteria are the current uncertainty estimates based on Dubovik et al. [2000]. 

701 
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703 Figure 1 . Distribution of the AERONET sites based on the dominant particle type. Sites were 

704 selected based on data volume, geographic location, and primary aerosol source region. Other 

705 dominant particle types (e.g., sea salt and biogenic aerosols) were not considered due to low 

706 aerosol loading conditions (T44o nm <0.4), which was a limiting threshold for AERONET Version 

707 2, Level 2.0 aerosol absorption retrievals [Dubovik et al., 2002; Holben et al., 2006]. 
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710 Figure 2. Spectral single scattering albedo averages were grouped by dominant aerosol particle 

711 category for T44o nm >0.4 using AERONET Version 2, Level 2.0 data. The plots utilize second 

712 order polynomial fit. 

713 
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Figure 3. Similar to Figure 2, except the 
spectral single scattering albedo averages for the 
Mixed category were grouped by fine mode 
fraction of AOD (r| 55 o n m) using the ranges 0.0- 
0.33 for coarse mode dominated particles (a), 
0.33-0.66 for mixed size particles (b), and 0.66- 
1.0 for fine mode dominated particles (c). 
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Figure 4. Absorption aerosol optical depth (x a bs) and absorption Angstrom exponent (a a b s ) 
averages were grouped by dominant aerosol particle category for x 44 o nm >0.4 using AERONET 
Version 2, Level 2.0 data. The plots use the power law fit and slopes of these lines are the a a b S 
(440-870 nm or 440-1020 nm) listed adjacent to the legend in each plot. 
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Figure 5. Similar to Figure 3, except i a b s and a a bs 
averages for the Mixed category were grouped by 
fine mode fraction of the AOD (r| 5 5 o n m) using 
ranges of 0.0-0.33 for coarse mode dominated 
particles (a), 0.33-0.66 for mixed size particles (b), 
and 0.66-1.0 for fine mode particles (c). 
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Figure 6. Absorption Angstrom exponent (a a b s ) frequency distribution for each dominant 
aerosol particle category using AERONET Version 2, Level 2.0 data. Approximately 10% of the 
a a bs retrievals (22% for Urban/Industrial) were below 1.0 or A,' 1 dependence. 
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743 Figure 7. Relative number density plots for the absorption Angstrom exponent (440-870 nm) 

o 

744 and extinction Angstrom exponent (440-870 nm) relationship based on dominant aerosol type 

745 using AERONET Version 2, Level 2.0 data. The color scale represents the relative density of 

746 points in each aerosol type partitioned data set, where orange to red colors (levels -45-64) 

747 indicate the highest number density based on the Voronoi tessellation. 
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749 Figure 8. Similar to Figure 7, except for the absorption Angstrom exponent (440-870 nm) and 

750 fine mode fraction of the aerosol optical depth (550 nm) relationship. 
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Figure 9. Similar to Figure 7, except for the single scattering albedo (440 nm) and the extinction 

o 

Angstrom exponent (440-870 nm) relationship. 
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755 Figure 10. Similar to Figure 7, except for the single scattering albedo (440 nm) and fine mode 

756 fraction of the aerosol optical depth (550 nm) relationship. 
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759 Figure 11. Weighted cluster averages were grouped for each aerosol type category and 

760 relationship using AERONET Version 2, Level 2.0 data. The Mixed category averages were 

761 calculated using a 0.8 extinction Angstrom exponent threshold between mainly small and mainly 

762 large particles. For the fine mode fraction of AOD, the Mixed category averages were calculated 

763 based on the 0.0-0.33, 0.33-0.66, and 0.66-1.0 ranges. For single scattering albedo plots, the 

764 Biomass Burning category was further partitioned by calculating averages using a single 

765 scattering albedo threshold of 0.90 to produce two sub-clusters (dashed ellipses) observed in 

766 Figure 9 and Figure 10. 
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